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1. Introduction 

The pioneering work [31] of Voiculescu connected free probability with random matrices, 
as one of the most prominent examples for a noncommutative probability space is the space 
of Hermitian N x N matrices. On one hand, the law of the sum of two free random variables 
with laws fia and is given by the free additive convolution ffl On the other 
hand, in case of Hermitian matrices, the law can be identihed with the distribution of the 
eigenvalues. Thus the free additive convolution computes the eigenvalue distribution of 
the sum of two free Hermitian matrices. However, freeness is characterized by an infinite 
collection of moment identities and cannot easily be verified in general. A fundamental direct 
mechanism to generate freeness is conjugation by random unitary matrices. More precisely, 
two large Hermitian random matrices are asymptotically free if the unitary transfer matrix 
between their eigenbases is Haar distributed. The most important example is when the 
spectra of the two matrices are deterministic and the unitary conjugation is the sole source of 
randomness. In other words, if A = and B = are two sequences of deterministic 
N X N Hermitian matrices and 17 is a Haar distributed unitary, then A and UBU* are 
asymptotically free in the large N limit and the eigenvalue distribution oi A + UBU* is 
given by the free additive convolution fj,A ffl /rs of the eigenvalue distributions of A and B. 

Since Voiculescu’s first proof, several alternative approaches have been developed, see e.g., 
[11, 16, 29, 30], but all of them were global in the sense that they describe the eigenvalue 
distribution in the weak limit, i.e., on the macroscopic scale, tested against iV-independent 
test functions (to fix the scaling, we assume that A^^'i and B^^'> are uniformly bounded). 


*Partially supported by ERC Advanced Grant RANMAT No. 338804. 
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The study of a local law, i.e., identification of the eigenvalue distribution of A+UBU* with 
the free additive convolution below the macroscopic scale, was initiated by Kargin. First, he 
reached the scale (log in [25] by using the Gromov-Milman concentration inequality 

for the Haar measure (a weaker concentration result was obtained earlier by Chatterjee [14]). 
Kargin later improved his result down to scale in the bulk of the spectrum [26] by 

analyzing the stability of the subordination equations more efficiently. This result was valid 
only away from finitely many points in the bulk spectrum and no effective control was 
given on this exceptional set. Recently in [1], we reduced the minimal scale to by 

establishing the optimal stability and by using a bootstrap procedure to successively localize 
the Gromov“Milman inequality from larger to smaller scales. Moreover, our result holds in 
the entire bulk spectrum. In fact, the key novelty in [1] was a new stability analysis in the 
entire bulk spectrum. 

The main result of the current paper is the local law for H = A + UBU* down to the 
scale for any 7 > 0. Note that the typical eigenvalue spacing is of order N~^, a 

scale where the eigenvalue density fluctuates and no local law holds. Thus our result holds 
down to the optimal scale. 

There are several motivations to establish such refinements of the macroscopic limit laws. 
First, such bounds are used as a priori estimates in the proofs of Wigner-Dyson-Mehta 
type universality results on local spectral statistics; see e.g., [21, 12, 20, 27] and references 
therein. Second, control on the diagonal resolvent matrix elements for some rj = Im z implies 
that the eigenvectors are delocalized on scale 77 “^; the optimal scale for rj yields complete 
delocalization of the eigenvectors. Third, the local law is ultimately related to an effective 
speed of convergence in Voiculescu’s theorem on the global scale [26, 1]. 

Our proof has three major ingredients. First, we use a partial randomness decomposition 
of the Haar measure that enables us to take partial expectations of Ga, the diagonal matrix 
elements of the resolvent G{z) = {H — z)~^ at spectral parameter z G C+. Exploiting 
concentration only for the partial randomness surpasses the more general but less flexible 
Gromov-Milman technique. Second, to compute the partial expectations of Gu, we establish 
a new system of self-consistent equations involving only two auxiliary quantities. Keeping in 
mind, as a close analogy, that freeness involves checking infinitely many moment conditions 
for monomials oi A, B and U, one may fear that an equation for G involves BG, whose 
equation involves BGB etc., i.e., one would end up with an infinite system of equations. 
Surprisingly this is not the case and monitoring two appropriately chosen quantities in 
tandem is sufficient to close the system. Third, to connect the partial expectation of Ga 
with the subordination functions from free probability, we rely on the optimal stability result 
for the subordination equations obtained in [ 1 ]. 

One prominent application of our work concerns the single ring theorem of Guionnet, 
Krishnapur and Zeitouni [22] on the eigenvalue distribution of matrices of the form UTV, 
where T is a fixed positive definite matrix and U, V are independent Haar distributed. Via 
the hermitization technique, the current proof of the local law for the addition of random 
matrices can also be used to prove a local version of the single ring theorem. This approach 
was demonstrated recently by Benaych-Georges [ 8 ], who proved a local single ring theorem 
on scale (logiV)“^/^ using Kargin’s local law on scale (logV)“^/^. The local law on the 
optimal scale N~^ and throughout the entire bulk spectrum is a key ingredient to prove the 
local single ring theorem on the optimal scale. The details are deferred to our separate work 
in preparation [ 2 ]. 

1.1. Notation. The following definition for high-probability estimates is suited for our pur¬ 
poses, which was first used in [18]. 
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Definition 1.1. Let 

X ={X<^^\v) : N €N,v€V^^'>), Y = (V’^^'>{v) : N €N,v €V‘^^'>) (1.1) 

be two families of nonnegative random variables where is a possibly N-dependent pa¬ 

rameter set. We say that Y stochastieally dominates X, uniformly in v, if for all (small) 
e > 0 and (large) D > 0, 

sup P fxW(u) > mY<^^'>{v)] < N-^ , (1.2) 

„ev(") V / 

for sufficiently large N > No{e,D). If Y stochastically dominates X, uniformly in v, we 
write X ^ Y. 

We further rely on the following notation. We use the symbols O(-) and o( •) for the 
standard big-0 and little-o notation. We use c and C to denote strictly positive constants 
that do not depend on N. Their values may change from line to line. For a, b > 0, we write 
Q ^ ^ ^ if there is C > 1 such that a < Cb, a > C~^b respectively. 

We use bold font for vectors in and denote the components as n = (ui, ..., uat) G . 
The canonical basis of is denoted by For v,w G C^, we write v*w for the scalar 

product We denote by ||n ||2 the Euclidean norm and by ||v||oo = max^ |ni| the 

uniform norm of n G C'^. 

We denote by Mm(C) the set oi N x N matrices over C. For A G MnIC), we denote 
by ll^ll its operator norm and by ||A ||2 its Hilbert-Schmidt norm. The matrix entries of 
A are denoted by Aij = e*Aej. We denote by tr^ the normalized trace of A^ i.e., ivA = 
X v,w G , the rank-one matrix vw* has elements {vw*)ij = {viWj). 

Let g = {gi,... ,gN) be a real or complex Gaussian vector. We write g ^ A/r(0, ct^In) if 
gi, ■ ■ ■, gN are independent and identically distributed (i.i.d.) 1V(0, cr^) normal variables; and 
we write g ^ A/c(0, u'^In) if gi, - ■ ■ ,gN are i.i.d. A^c(0, cr^) variables, where gi ^ Nc{0, cr'^) 
means that Reg^ and Iragi are independent 1V(0, ^) normal variables. 

Finally, we use double brackets to denote index sets, i.e., 

|ni,n2l := [ni,n2] nZ, 

for ni, 712 G M. 


2. Main results 

2.1. Free additive convolution. In this subsection, we recall the definition of the free 
additive convolution. This is a shortened version of Section 2.1 of [1] added for completeness. 

Given a probability measure* ^ on K its Stieltjes transform, m^, on the complex upper 
half-plane C"*" := {z G C : Imz > 0} is defined by 

.»„«:= ( 2 . 1 ) 

X — z 

Note that rui^ : C+ —>■ C"*" is an analytic function such that 

lim ir]mfj,(ig) = —1. (2.2) 

r]y^OO 

Gonversely, if m : C"*" —>■ C“*" is an analytic function such that lim^/to^ ir]m{ir]) = 1, then m 
is the Stieltjes transform of a probability measure g, i.e., m{z) = m^{z), for all 2 G C+. 


*A11 probability measures considered will be assumed to be Borel. 
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We denote by the negative reciprocal Stieltjes transform of i.e., 




1 

m^(z) ’ 


zeC+. 


(2.3) 


Observe that 


lim 

r]/'oo ITJ 


= 1 , 


(2.4) 


as follows from (2.2). Note, moreover, that Ff^ is an analytic function on C"*" with nonnegative 
imaginary part. 

The free additive convolution is the symmetric binary operation on probability measures 
on K characterized by the following result. 


Proposition 2.1 (Theorem 4.1 in [6], Theorem 2.1 in [15]). Given two probability measures, 
p,i and p, 2 , on K, there exist unique analytic functions, wi ,uj 2 '■ C"*" —>■ C"*", such that, 

(i) for all z € C’*', Ima;i(z), Ima; 2 ( 2 ;) > Im 2 , and 


lim 

rj/^oo 


i-q 


(ii) for all z G C’*', 


lim 

rj/^oo 


c^2(i??) 

iq 


= 1 ; 


UJl{z) +UJ2{z) - Z = Ff^^{uJ2{z)) . 


(2.5) 

( 2 . 6 ) 


It follows from (2.5) that the analytic function F : C’*' —>■ C’*' defined by 


F{z) := F^,{uj2{z)) = Ff,^{uJi{z)), 


(2.7) 


satisfies the analogue of (2.4). Thus F is the negative reciprocal Stieltjes transform of a 
probability measure /r, called the free additive convolution of q,i and /j, 2 , usually denoted by 
/r = /ri ffl /i 2 . The functions uii and uj 2 of Proposition 2.1 are called subordination functions 
and m is said to be subordinated to , respectively to . Moreover, observe that uji 
and UJ 2 are analytic functions on with nonnegative imaginary parts. Hence they admit 
the Nevanlinna representations 

f 1 ~\~ zx 

uij{z) = aujj+z+ -d£i^^.(a;), j = l,2, z S C+, (2.8) 

JR X — Z 

where € R and are finite Borel measures on K. For further details and historical 
remarks on the free additive convolution we refer to, e.g., [32, 23]. 

Choosing qi as a single point mass at 6 € R and fi 2 arbitrary, it is straightforward to check 
that fii ffl ^2 is ^*2 shifted by b. We exclude this uninteresting case by assuming hereafter 
that /ii and ^2 are both supported at more than one point. For general pLi and ^ 2 , the 
atoms of pLi ffl ^2 are identified as follows. A point c G R is an atom of /ii ffl ^ 2 , if and only if 
there exist a, 6 € R such that c = a + b and /ii({a}) + ti 2 {{b}) > 1; see [Theorem 7.4, [10]]. 
Properties of the continuous part of pLi ffl fi 2 may be inferred from the boundary behavior 
of the functions oji and UJ 2 - For simplicity, we restrict the discussion to compactly 

supported probability measures. 


Proposition 2.2 (Theorem 2.3 in [3], Theorem 3.3 in [4]). Let fii and ^2 be compactly 
supported probability measures on R none of them being a single point mass. Then the 
functions uii, UJ 2 ■ C'*' —?► C"*" extend continuously to R. 

Belinschi further showed in Theorem 4.1 in [4] that the singular continuous part of /ii ffl ^2 
is always zero and that the absolutely continuous part, (/ii ffl M 2 )'^'^i of ffl fi 2 is always 
nonzero. We denote the density function of {fii ffl ^ 2 ^^ by 
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We are now all set to introduce our notion of regular hulk, Mi Informally, 

we let be the open set on which ffl /j ,2 has a continuous density that is strictly 

positive and bounded from above. For a formal definition we first introduce the set 

•= int I supp {gLi ffl ^ 2 )’"'" \ {x G R : lim + i??) = 0} I. (2.9) 

f »?\.o J 

Note that does not contain any atoms of fj,i ffl /i 2 . By the Luzin-Privalov theorem 

the set {a: G M : lim^j^^o ( 2 ; + i^) = 0} has Lebesgue measure zero. In fact, an 

even stronger statement applies for the case at hand. Belinschi [5] showed that if a: G M 
is such that lim^\^o ^/xiffl/i 2 ( 2 ^ + = 0 i then it must be of the form x = a + b with 

Mi({®}) + M 2 ({^}) > I, a, & G R. There could only be finitely many such x, thus 
must contain an open non-empty interval. 

Proposition 2.3 (Theorem 3.3 in [4]). Let fii and p ,2 be as above and fix any x G 
Then uii, 012 : C"*" —?► C'*' extend analytically around x. Thus the density function 

/^iH ^2 analytic in wherever positive. 

The regular bulk is obtained from by removing the zeros of inside 

Definition 2.4. The regular bulk of the measure yii ffl /r 2 is the set 

^/ilffl/i2 ■ ^/ilffl/i2 \ {2" S ^/iiffl/i2 ■ /aIiEBaI2 ( 2 ^) ■ (2-10) 

Note that is an open nonempty set on which fii ffl ^2 admits the density 

The density is strictly positive and thus (by Proposition 2.3) real analytic on 

2.2. Definition of the model and assumptions. Let A = and B = be 

two sequences of deterministic real diagonal matrices in Mjv(C), whose empirical spectral 
distributions are denoted by pi a and piB, respectively. More precisely, 

AT AT 

TA-=J^'^Sai, bl'B ■= , (2.II) 

i—1 i—1 

with A = diag(ai), B = diag( 6 i). The matrices A and B actually depend on TV, but we omit 
this from our notation. Proposition 2.1 asserts the existence of unique analytic functions oja 
and ujb satisfying the analogue of (2.5) such that, for all z G C+, 

Ff,^{ujB{z)) = Ffj,^{u}A{z)) , u}a{z)+ujb{z) - z = Fij,^{ujb{z)) . ( 2 . 12 ) 

We will assume that there are deterministic probability measures pLa and on R, neither 
of them being a single point mass, such that the empirical spectral distributions p-A and pB 
converge weakly, as iV —>■ 00 , to pa and pp, respectively. More precisely, we assume that 

dL(MA,Mc<)+ dL(/iSiM/3)0, (2.13) 

as TV —>• 00 , where du denotes the Levy distance. Proposition 2.1 asserts that there are 
unique analytic functions Wq, top satisfying the analogue of (2.5) such that, for all z G C’*', 

= Pufli^^aiz)) , UJa{z) + U}p{z) - Z = Ffj,^{uJp{z)) . (2.14) 

Proposition 4.13 of [9] states that dL(MA ffl PB,Pa ffl Pp) < dL(MAjMa) + ‘ih{pB, pp), Le., 
the free additive convolution is continuous with respect to weak convergence of measures. 

Denote by U (TV) the unitary group of degree TV. Let U €U (TV) be distributed according 
to the Haar measure (in short a Haar unitary), and consider the random matrix 

H = :=A + UBU* . 


(2.15) 
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Our results also holds for the real case when U is Haar distributed on the orthogonal group, 
0{N), of degree N. Throughout the main part of the paper the discussion will focus on the 
unitary case while the orthogonal case is addressed in Appendix A. 

We introduce the Green function, Gh, of H and its normalized trace, mn, by 

Gh{z ):=-^—, := trG'//(z), z G C~^ . (2.16) 

H — z 

For simplicity, we frequently use the notation G{z) instead of Gh{z) and we write Gij{z) = 
{GH)ij{z) for the (*, j)th matrix element of G{z). 


2.3. Main results. For a,b > 0, b > a, and I C K, let 

iSi(a, b) := {z = E + irj £ C"*" : EGl,a<rj<b}, (2-17) 

In addition, for brevity, we set, for any given 7 > 0, 

Vm = Vmil) ■= (2.18) 

The main results of this paper are as follows. 


Theorem 2.5. Let o,nd up be two compactly supported probability measures on R, and 
assume that neither is supported at a single point and that at least one of them is supported 
at more than two points. Assume that the sequences of matrices A and B in (2.15) are such 
that their empirical eigenvalue distributions pA and ps satisfy (2.13). Let I C be 

a nonempty compact interval. 

Then, for any fixed 7 > 0, the estimates 


max 

l<i<Af 


G,fz) 


max 


G^jiz) 


1 


Qi - UJb{z) 


1 


1 


(2.19) 

( 2 . 20 ) 


and 


m//(2;)-m^^EB^3(z) 






( 2 . 21 ) 


hold uniformly on Sx(r]m, 1 ) (see where rj = Imz and rj^n is given in (2.18). 


Remark 2.1. The assumption that neither of pa and pp is a point mass, ensures that the 
free additive convolution is not a simple translate. The additional assumption that at least 
one of them is supported at more than two points is for the brevity of the exposition here. 
In Appendix B, we present the corresponding result for the special case when pa and pp are 
both convex combinations of two point masses. 


Remark 2.2. We recall from Lemma 5.1 and Theorem 2.7 of [1] that, under the conditions 
of Theorem 2.5, there is a finite constant G such that 

max max I \ujAiz) - uJa{z)\, \ujb{z) - ujp{z)\, |m^^ffl^3(z) - (z)| I 

.2 G»Sz(0 ,1) ^ ^ 

< G {dB{pA-, Pa) + dxipB, pp)) , (2.22) 

i.e., the Levy distances of the empirical eigenvalue distributions of A and B from their 
limiting distributions control uniformly the deviations of the corresponding subordination 
functions and Stieltjes transforms. Note moreover that maXj,g 5 ^(Q 1 ) |w^^h^^( 2 ;)| < 00 by 
compactness of X and analyticity of Thus the Stieltjes-Perron inversion formula 

directly implies that {pA ffl Pb)^"^ has a density, inside X and that 

^^^\ftJ.Aap.Bix) - fucaij.f3{x)\ < G{dx{pA,Pa) + dB{pB,pp)) , 


(2.23) 
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for N sufficiently large. In particular, since I C we have I C for N 

sufficiently large, i.e., we use, for large TV, /ia to locate (an interval of) the regular bulk 
of ffl Mb- Finally, combining (2.22) and (2.19), we get 


max 

l<i<N 


Gu{z) 


ai-ujp{z) y^Nrj 




+ dL(M7l,Ma) + dL(MB,M/3) : 


uniformly on 1), where rj = Imz. Averaging over the index i, we get the correspond¬ 
ing statement for \mH — with the same error bound. 


Remark 2.3. Note that assumption (2.13) does not exclude that the matrix H has outliers 
in the large N limit. In fact, the model H = A + UBU* shows a rich phenomenology when, 
say, A has a finite number of large spikes; we refer to the recent works in [7, 13, 26]. 


Let Ai,..., Aat be the eigenvalues of H, and iti,..., mtv be the corresponding £^-normalized 
eigenvectors. The following result shows complete delocalization of the bulk eigenvectors. 


Theorem 2.6 (Delocalization of eigenvectors). Under the assumptions of Theorem 2.5 the 
following holds. LetX C be a compact nonempty interval. Then 


max 


uA 


-< 


1 

7 ^ ■ 


(2.24) 


2.4. Strategy of proof. In this subsection, we informally outline the strategy of our proofs. 
Throughout the paper, without loss of generality, we assume 


trA = trB = 0. 


(2.25) 


For brevity, we use the shorthand for the Stieltjes transform of ma H Ms- 

We consider first the unitary setting. Let 

H:=A + UBU*, n-.= U*AU + B, (2.26) 


and denote their Green functions by 

G(z) = (i7-z)-\ g{z) = {n-z)-\ zec+. (2.27) 


We write z = E + iri G C'*', with E gM. and p > 0, for the spectral parameter. In the sequel 
we often omit z € C"*" from the notation if no confusion can arise. Recalling (2.16), we have 

mniz) = trG(z) = trt/(z), z G C*". 

For brevity, we set 

A:=U*AU, B:=UBU*. 


The following functions will play a key role in our proof. 


Definition 2.7 (Approximate subordination functions). 


(jj\{z) :=z- 


tr Ag{z) 
mniz) 


uj%{z) := z - 


trRG(z) 

mniz) 


z G C+ . (2.28) 


Notice that the role of A and B are not symmetric in these notations. By cyclicity of the 
trace, we may write 

ujI(z) = z - , ZGC+. (2.29) 

mH[z) 

We remark that the approximate subordination functions defined above are slightly different 
from the candidate subordination functions introduced in [26] which were later used in [1]. 
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The functions oj\[z) and w^(z) turn out to be good approximations to the subordination 
functions uja{z) and ujb{z) of (2.12). A direct consequence of the definition in (2.28) is that 

—= z - 0J%{z) - Uj%{z), zec+. (2.30) 

mH[z) 

Having set the notation, our main task is to show that 

Gii(z) = (tti — uj^iz)) +0^^—^==), z € 5x(?7m, 1) j (2-31) 

where we focus, for simplicity, on the diagonal Green function entries only. 

We first heuristically explain how (2.31) leads to our main result in (2.19). A key input is 
the local stability of the system (2.12) established in [1]; see Subsection 3.3 for a summary. 
Averaging over i in (2.31), we get 

mH{z) = mA{uJB(.z))+ 0^(^-^y (2.32) 

Replacing H by H, we analogously get 

niHiz) = mB{uj‘^{z))+ 0^(^-^y (2.33) 

according to (2.29). Substituting (2.30) into (2.32) and (2.33) we obtain the system 

— 2 + 0^ (yW?)) ’ 

which is a perturbation of (2.12). Using the local stability of the system (2.12), we obtain 

\^a{z) - t^A{z)\ , \u}%{z) - u}b{z)\ . (2.34) 

Plugging this estimate back into (2.31) we get (2.19). The full proof of this step is accom¬ 
plished in Section 7. 

We now return to (2.31). Its proof relies on the following decomposition of Haar measure 
on the unitary group given, e.g., in [17, 28]. For any fixed i £ |1, A^J, any Haar unitary U 
can be written as 

U = (2.35) 

Here Ri is the Householder reflection (up to a sign) sending the vector to Vi, where 
Vi £ is a random vector distributed uniformly on the complex unit {N — l)-sphere, and 
6i £ [0, 27r) is the argument of the ith coordinate of Vi. The unitary matrix has as its 
zth column and its (f,i)-matrix minor (obtained by removing the fth column and ith row) 
is Haar distributed on [/(N — 1); see Section 4 for more detail. 

The gist of the decomposition in (2.35) is that the Householder reflection Ri and the 
unitary are independent, for each fixed i £ |1, A^]. Hence, the decomposition in (2.35) 
allows one to split off the partial randomness of the vector Vi from U. 

The proof of (2.31) is divided into two parts: 

{i) Concentration of Ga around the partial expectation i.e., 

\Gu-^.AG,i]\ <^=. 

y/Nr] 
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[ii) Computation of the partial expectation i.e., 

[Gii{z)] - (a, - uj%{z))~^\ -< . 

To prove part (i), we resolve dependences by expansion and use concentration estimates 
for the vector Vi. This part is accomplished in Section 5. 

Part {ii) is carried out in Section 6. We start from the Green function identity 

(a* - z)Gu{z) = -{BG{z))u + 1 ■ (2.36) 


Taking the expectation of (2.36) and recalling the definition of the approximate subor¬ 
dination function uj'^{z) in (2.28), it suffices to show that 

to prove (2.31). Denoting := and setting, for z G C+, 

Sf(z) := e‘®“<il<*>G(z)ei, Tf(z) := e'^^v*G(z)ei, 


we will prove that 

E„, [(5g(z))..] = -E„. [5»(z)] + (^) . (2.37) 

Hence, it suffices to estimate E„. [S'?] instead. Approximating by a Gaussian vector 

and using integration by parts for Gaussian random variables, we get the pair of equations 

E„. [S“] = tr {BG) (E„. [sf] - h E„, [rf]) + tr (HGH) {Gu + E„, [xf]) + (^) > 

E„, [Tf] = tr (G) (E., [Sf] - h E„, [Tf]) + tr (BG) (G.. + E., [T#]) + (^) , 


where we dropped the z-argument for the sake of brevity; see (6.23) and (6.24) for precise 
statements with slightly modified Sj^ andTj'^. Solving E„. [S|*] from the above two equations, 
we arrive at 

V 2 


E.. [sf] = - 


tr (BG) 


trG 


G,, 


tr(HG) - (tr (HG))" 


tr G 


tr {BGB) 


(G,,+E„, [Tf]) 


O. 




(2.38) 


Returning to (2.37), we also obtain, using concentration estimates for {BG)ii (which follow 
from the concentration estimates of Gu established in part (i) and (2.36)), that 

N 


l^E„.[S#]+trHG 




1 




(2.39) 


Thus, averaging (2.38) over the index i and comparing with (2.39), we conclude that 


tr {BG) - (tr (HG))' 


trG 


-h tr (BGB) 


-< 


y^Nrj 


Plugging this last estimate back into (2.38), we eventually find that 
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which together with (2.37) and (2.36) gives us part (ii). This completes the sketch of the 
proof for the unitary case. The proof of the orthogonal case is similar. The necessary 
modifications are given in Appendix A. 


3. Preliminaries 

In this section, we first collect some basic tools used later on and then summarize results 
of [1]. In particular, we discuss, under the assumptions of Theorem 2.5, stability properties 
of the system (2.12) and state essential properties of the subordination functions uja and ujb- 


3.1. Stochastic domination and large deviation properties. Recall the definition of 
stochastic domination in Definition 1.1. The relation ^ is a partial ordering: it is transitive 
and it satisfies the arithmetic rules of an order relation, e.g., if Xi -< Yi and X 2 -< I 2 then 
Xi + X 2 <Yi + Y 2 and X 1 X 2 -< Y 1 Y 2 . Further assume that $(u) > N~^ is deterministic 
and that Y{v) is a nonnegative random variable satisfying E[F(i;)]^ < N'^ for all v. Then 
Y{v) -< d>(i;), uniformly in u, implies E[y(u)] ^ d>(t;), uniformly in v. 

Gaussian vectors have well-known large deviation properties. We will use them in the 
following form whose proof is standard. 

Lemma 3.1. Let X = {xij) S Mv(C) he a deterministic matrix and let y = (yi) £ be a 
deterministic complex vector. For a Gaussian random vector g = {gi ,..., gjq) £ A/r(0, 
or Afc{0,o-^lN), we have 

\y*g\ F a\\y\\ 2 , \g*Xg - a'^NtrXl ^ a‘^\\X\\ 2 . (3.1) 


3.2. Rank-one perturbation formula. At various places, we use the following funda¬ 
mental perturbation formula: for a,/3 £ and an invertible D £ MAr(C), we have 


{D + af3*) ^ = D-i 


D-^al3*D-^ 

1 -f/3*D-ia ’ 


(3.2) 


as can be checked readily. A standard application of (3.2) is recorded in the following lemma. 


Lemma 3.2. Let D £ Mjv(C) he Hermitian and let Q £ he arbitrary. Then, for 

any finite-rank Hermitian matrix R £ Mn{C), we have 


ti{Q[D + R-z) ^^-tT[Q{D-z) 


< 


rank(i?)||(5|| 


z = E + iyGC+. (3.3) 


= ± 


Proof. Let z £ C"*" and a. £ . Then from (3.2) we have 

1 a.*{D — z)~^Q{D — z)~^a. 
N 1 F a*{D — z)~^a. 

Iigil \\{D-z)-^a\\l 


tr ^g(Zl ± aa* — z) — tr ^g(D — z) 


We can thus estimate 


tr 


(q( 


D ± OLOL* - 


z) -tv{Q{D-z) 


< 


< 


N 11 ± q:*(D — z)“^q:| 
IIQII a.*lm{D — z)~^OL 
Nrj 11 ± q:*(D — z)“^q:| 

IIQII 


Nrj 


(3.4) 


(3.5) 


Since R = R* € M]\j{<C) has finite rank, we can write i? as a finite sum of rank-one Hermitian 
matrices of the form iaa*. Thus iterating (3.5) we get (3.3). □ 
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3.3. Local stability of the system (2.12). We first consider (2.6) in a general setting: 
For generic probability measures /ri,/r 2 , let be given by 




-OJ1-UJ2+ Z \ 
-UJ1-UJ2+ Z ) ’ 


(3.6) 


where are the negative reciprocal Stieltjes transforms of ^ 1 , [i 2 \ see (2.3). Consid- 

ering ^i ,//2 as fixed, the equation 


i/J'2 (^17 ^2 ; ■2^) O7 


(3.7) 


is equivalent to (2.6) and, by Proposition 2.1, there are unique analytic functions Wi,aj 2 : 
C+ —>■ C+, z !->• a;i(z), u; 2 (z) satisfying (2.5) that solve (3.7) in terms of z. Choosing /ri = 

= /.t /3 equation (3.7) is equivalent to (2.14); choosing fii = ^2 = Mb it is equivalent 

to (2.12). When no confusion can arise, we simply write 4> for (wi, W 2 , z). 

We call the system (3.7) linearly S-stable at (wi,W 2 ) if 


rMi.A‘2 (wi, W 2 ) 


-1 F'^^{u;2) 

f;^(c.i)-i -1 



(3.8) 


for some positive constant S. In particular, the partial Jacobian matrix of (3.6) given by 


D$(a;i,W 2 ) := 


94* 9$ 

— (cui,u;2,z),—(a;i,a.2,.)) = 


-1 

F'M)-^ 


-1 


has a bounded inverse at {u!i,ui 2 )- Note that D$(a;i,W 2 ) is independent of z. The implicit 
function theorem reveals that, if (3.7) is linearly S'-stable at (a;i,W 2 ), then 


max (z)| < 2S , 

z 6 Sx( 0 ,l) 


max \ujUz)\<2S. 

26 Si( 0 , 1 ) 


(3.9) 


In particular, uii and W 2 are Lipschitz continuous with constant 25'. A more detailed analysis 
yields the following local stability result of the system (wi, a; 2 , 2 ) = 0. 


Lemma 3.3 (Proposition 4.1, [1]). Fix Zq € C“'". Assume that the functions oji, 102 , ri, 
r 2 ■ C"*" —>■ C satisfy Imwi(zo) > 0, lmw 2 (zo) > 0 and 


(a;i(zo),W 2 (zo),Zo) = r(zo), 


(3.10) 


where r(z) '■= (ri(z),r 2 (z))^. Assume moreover that there is S G [0,1] such that 

|wi(zo) - a;i(zo)| < J, |w 2 (zo) - a; 2 (zo)| < <5, (3.11) 


where uji{z), ui 2 {z) solve the unperturbed system (wi,a; 2 , z) = 0 with lma;i(z) > Imz 

and Ima; 2 (z) > z, z G C"*". Assume that there is a constant S such that is linearly S-stable 
at (a;i(zo),a; 2 ( 2 :o)); assume in addition that there are strictly positive constants K and 
k with k > 6 and with k^ > SKS such that 


Then 


k <1 tiiijJi{zq) < K, fc < lma;2(zo) < AT. 


(3.12) 


|a;i(zo) - wi(zo)| < 25||r(zo)||2 , |w 2 (zo) - a; 2 (zo)| < 25||r(zo)||2 ■ (3.13) 

In Section 7, we will apply Lemma 3.3 with the choices mi = Ta and pL 2 = ps- We 
thus next show that the system (w^i, w_b, z) = 0 is 5-stable, for all z G 5x(0,1), and 

that (3.12) holds uniformly on 5i(0,1); see (2.17) for the definition. 
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Lemma 3.4 (Lemma 5.1 and Corollary 5.2 of [1]). Let ha, Hb be the probability measures 
from (2.11) satisfying the assumptions of Theorem 2.5. Let uja,^b denote the associated 
subordination functions of (2.12). Let T be the interval in Theorem 2.5. Then for N suf¬ 
ficiently large, the system (w^, ws,z) = 0 zs S-stable with some positive constant S, 

uniformly on 5x(0,1). Moreover, there exist two strictly positive constants K and k, such 
that for N sufficiently large, we have 


max |a;^(z)| < K , 
^G5i(0.1) 

max \u}b(z)\ < K , 

^G5r(0.1) 

(3.14) 

min Imw^(z) > k, 

2 G 5 x( 0 , 1 ) 

min ImuJBiz) > k. 
zeSi{Q,i) 

(3.15) 

Remark 3.1. Under the assumptions of Lemma 3.4, the estimates in (3.15) can be extended 
as follows. There is fc > 0 such that 

min (lmw4(z) — Imz) > , 

^G 5 x( 0 , 1 ) 

min (Ima;B(z) — Imz) > fc . 
^eSz(o,i) 

(3.16) 


This follows by combining (3.15) with the Nevanlinna representations in (2.8). 

We conclude this section by mentioning that the general perturbation result in Lemma 3.3 
combined with Lemma 3.4, can be used to prove (2.22). We refer to [1] for details. 


4. Partial randomness decomposition 


We use a decomposition of Haar measure on the unitary groups obtained in [17] (see 
also [28]): For a Haar distributed unitary matrix U = Un, there exist a random vector 
Vi = {vii,... ,vin), uniformly distributed on the complex unit {N — l)-sphere := 

{x G : x*x = 1}, and a Haar distributed unitary matrix G U{N — 1), 

which is independent of Ui, such that one has the decomposition 


U = -e‘®i 




where 


ri 


:= V2-, 


ei + e 

ei + e-‘®iuil]2 


R, := 1 — r,r 


I'l 


(4.1) 


and where 0i is the argument of the first coordinate of the vector Vi. More generally, for 
any i G |1,A^], there exists an independent pair with Vi a uniformly distributed 

unit vector Vi and with G U{N — 1) a Haar unitary, such that one has the decomposition 


U = , 


:= V 2 , 


Gi + e '^'Vi 


R, —L-r.r* 


(4.2) 


where 1/^*^ is the unitary matrix with as its zth column and C/® as its (z,z)-matrix minor. 
With the above notation, we can write 

H = A + Ri&^R* , 


for any i G |1, A^], where we introduced the shorthand notation 

We further define 

ijd) := H (z) := (i7<*> - z)-\ 

Note that and are independent of Vi. 


(4.3) 


zGC+. 


(4.4) 
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It is well known that for a uniformly distributed unit vector Vi € C^, there exists a 
Gaussian vector = (gn,- ■ ■ , giN) ^ A/c(0, N~^I) such that 


V. = 


\\a^\\2 


By definition, 9i is also the argument of gu. Set 


9ik — 6 Qiiz 


k ^ i, 


(4.5) 


and introduce an iVc(0, N~^) variable gu which is independent of the unitary matrix U and 
oiQu Then, we denote g, := {gu,..., gi^) and note ~ A/c(0, In addition, by 

definition, we have 


\gii\ gi' 


—Wi li/Ml ya I / 1 


In subsequent estimates for Gij, it is convenient to approximate by 

:= Bi + Qj 


(4.6) 


in the decomposition U = without changing the randomness of To esti¬ 

mate the precision of this approximation, we require more notation: Let 

Wi=W* := I - w,w* , b''^ = WiB^^W^. (4.7) 

Correspondingly, we denote 

A + G‘'^{z):={H‘'^ - z)-^, z€C+. (4.8) 


The following lemma shows that can be replaced by Wi in Green function entries at 
the expense of an error that is below the precision we are interested in. 


Lemma 4.1. Fix z = E + \g £ C"*" and choose indices i,j, k £ |1, iVJ|. Suppose that 



max||Gfcfc(z)|, |G,^*^( 2 :)|| ^ 1, 
max||g*GW(2)ej|,|g*i?WGW(z)ej|| ^ 1, 

(4.9) 

hold. Then 

\Gkj{z) G^^{z)\^^ 

(4.10) 

holds, too. 




Proof of Lemma 4-L Fix i,j,k £ [1, A^]. We first note that 

ri = w^ + SuBi -I- S2tg^ , 

where 


6ii 

52i '■ = 


V 2 


\ei + e-^^'Vih 


V2 


-1 + 


1 


\gii\ gi\ 

\B^ + e-'^^'Vi\\2 Wgih 


- 1 . 


\B^ + e-^^’Vi\\2 \\g,\\2 

By the strong concentration of the norms in (4.11) and gu^'gu Nc{Q, N~^), we have 


(4.11) 




1 

’ 




y/N ' 


(4.12) 
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Denote 


Ai := w,w* - Vir* . 

Fix now z € C"*". Dropping z from the notation, a first order Neumann expansion of the 
resolvent yields 

Gkj = eg - (G(A,5«IF, + A, + A,i3« A,)GW)^^. . (4.13) 

Observe that the second term on the right side of (4.13) is a polynomial in the terms 



gtG^^e, g*B«G«e,, 

e*5«G(bej-, 


Gki , 

elGgi, elGB^^'>g,, 

elGB^^'^e,, 

(4.14) 






with coefficients of the form , for some nonnegative integers fci, k 2 such that fci +^2 > 1. 

By assumption (4.9) and the fact = biSi, we further observe that the first four terms 

in (4.14) are stochastically dominated by one. The last four terms are also stochastically 
dominated by one as follows from the trivial fact e*B^'^^ei = bi and Lemma 3.1. The terms 
in the second line of (4.14) are stochastically dominated by 

|e*,GQ«x,| ^ ||Q«||||Gefe||2 < ^(G^G)^^ = -< ^ , (4.15) 

with = / or B^'‘\ and with Xi = Bi or g^, where the last step follows from (4.9). Note 
that the terms in the second line of (4.14) appear only linearly in (4.13). Hence, (4.12), (4.15) 
and the order one bound for the first and last four terms in (4.14) lead to (4.10). □ 

5. Concentration with respect to the vector 

In this section, we show that concentrates around the partial expectation EgJG^*^], 
where Eg. [ • ] is the expectation with respect to the collection (Re gtj , Im gij. Besides the 
diagonal Green function entries = e*G^®^ei the following combinations are of importance 

T,(z) :=g:GW(z)e,, 5,(z) := g^H^GWe,, z G C+. (5.1) 

The estimation of EgJG-*^], carried out in the Sections 6 and 7, involves the quantities Ti 
and Si- From a technical point of view, it is convenient to be able to go back and forth 
between T^, Si and their expectations Eg. [T^], Eg. Thus after establishing concentration 
estimates for G)/ in Lemma 5.1 below, we establish in Corollary 5.2 concentration estimates 
for Ti and Si where we also give a rough bounds on Ti, Si and related quantities. We need 
some more notation: for a general random variable X we define 

lEg.A := A-EgW. (5.2) 

The main task in this section is to prove the following lemma. 

Lemma 5.1. Suppose that the assumptions of Theorem 2.5 are satisfied and let ^ > 0. Fix 
z = E + iri € Sx{r]ni, 1) and assume that 

\G^i{z) - ^ N~i , \Gli’{z)-{ai-u}Biz))~^\^ N~i , (5.3) 
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uniformly in i G |1,A^]. Then 


max 




^/Nri ' 


(5.4) 


Proof of Lemma 5.1. In this proof we fix z G 5i(?7ni,l). Recall the definition of G^'^\z) 
in (4.4) and note that G^'^\z) is independent of Vi (or gf). It is therefore natural to expand 
(z) around G^®^ (z) and to use the independence between G^®^ (z) and in order to verify 
the concentration estimates. However, by construction, we have 

= TTTT-T’ (^-5) 

(^i ~T Oi — Z 

which may be as large as l/? 7 , depending on Uj, bi and z. To circumvent problems coming 
from instabilities in G^f (z), we may use a “regularization” trick to enhance stability in the 
ei-direction: instead of considering the Green function of = A + H*^®^ directly, we first 
consider the (z-dependent) matrix 

i7^®^(z) := A + H^®) - {bi + ujb{z) - z)eie* , (5.6) 

and define G{®>(z) := (H{®>(z) - z)-^ Note that ij{®>(z) is not symmetric, yet since 
Ima;B(z) > Imz by Proposition 2.1, G^®^(z) is in fact well-defined on the whole upper-half 
plane. Fix any j G |I,iV|. Using the rank-one perturbation formula (3.2), we get 


g(®)(^) = _ (f’^+^B(z)-z)Gj;\z)G};^(z) 

^ + {b^+^B{z)- z)G\^{z) 
Some algebra then reveals that 

:{*}/ 




Gg^(^) 


1 + (bi +u}b{z) - z)Gji\z) 


^ - ih+u^Bjz) - z)inEgJGW(z)] G);^(z) 

1 + + UJb {z) - z)Eg. [GI-^ (z)] 




By assumption (5.3) and identity (5.7), we have 


G^iz)- 


a* - ujb(z) 


-< N-i 


G\f^{z)- 


I 


^ N-i . 


(5.7) 


(5.8) 


(5.9) 


Qi-bi- 2ujb{z) + z 

Note that |Im(ai —6 ^ —2wb(z)-|-z)| > ImwB(z) > 0, thus both denominators are well sepa¬ 
rated away from 0 by Lemma 3.4, in particular G\f (z) and g|*^ (z) are uniformly bounded. 
We will prove below that 

(5.10) 

^/Nfj ^ ^ 

Setting j = i in (5.8) and expressing the denominator on the right side by using (5.9)-(5.10), 
we get 

Oi - ujb{z) 


5^JG«(z)]h 


I -f {b,+ujB{z) - z)Eg^[G\-\z)] - 


^ N-z . 


(5.II) 


Oi-bi- 2 ujb{z) + z 

In particular, together with Lemma 3.4 and Ima;B(z) > Imz, this implies that the absolute 
value of the denominator on the right side of (5.8) is bounded from below by some strictly 
positive constant. Thus, applying lEg. on both sides of (5.10), we obtain the concentration 
estimate in (5.4). 

In the rest of the proof, we verify (5.10). Consider next the matrix 
LfW(z) := H + H<®> - (6, + ujb{z) - z)e,e* , 


(5.12) 
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and let G[*1 (z) := — z)~^. Note that depends on z and ujb{z) and is thus not 

symmetric, yet since ln\u)B{z) > Imz, GW(z) is, similar to G^®^, well-defined on the whole 
upper-half plane. Note that 




1 

Qi - COb{z) 


< 1 


(5.13) 


since Imws is uniformly bounded from below on 1) by Lemma 3.4. 

We now expand G^®^(z) around G^'^\z) and use the independence among G'^'^'^{z) and 
For simplicity, we hereafter drop the z-dependence from the notation. We start with noticing 
that 


^{i} _ ^[4 = -WiW*B^^'> - + w,w*B^^'>WiW* 


= —WjW 

= W,s 




WiW, 




where we introduced 

Si := -B^'^'>w^ , ti := -( 5 ^®^ - w*B'''‘'>Wi l)w^ 
Iterating the rank-one perturbation formula (3.2) once, we obtain 


G{®> = XW - 


X[®lGm*X[®l 


l + w*X[H, ’ 

Taking the (j,j)th matrix entry in (5.16), we have 


X[®1 := gW - 


1 -b s*G[®l 


Wi 


■^{^} _ /^[^] 




— ^ 0 ^ I 

^3 




where we introduced 


(5.14) 

(5.15) 


(5.16) 


(5.17) 


and 


E, := {s*G^^w,) + {w*G^%) + {s*G^^w,){w*G^%) - {w*G^^w,){s*G^%), (5.18) 

:= -(e*GWt,)((«^:G[®le,) + {s*G^^w,){w*G^^e,) - {w*G^^w,)is*G^^e,)) 

- {e*G^^w,){{s*G^’^ej) + iw*G^%){s*G^^ej) - {s*G^^U){w*G^^ej)) . (5.19) 


We now rewrite (5.17) as 




- 


(Gl 


Id 


1 [’^i] 



(5.20) 


Since |gJ*^| ^ 1 (c./., (5.9)) and |g[*^| ^ 1 (c./., (5.13)), it suffices to verify the following 
statements to show (5.10): 

(i): |]nEgjS,]| ^ (ii): _^-^^l, (iii): | ^ ^. (5.21) 

We first show claim (f). Substituting the definitions in (5.15) into (5.18), we have 
S, = - w*B<^'>G^^Wi - w*G^’^B<^^w, + w*B<^'>G^^Wi w*G^’^B<^'>Wi 

- w*B<^^G^’^B<^'>w,+w*B<^^Wi w*G^’^w, 

- w*W'^w, w*B^^G^^w, - w*W'^w, w*G^^w, w*B^^'>G^^Wi. 


(5.22) 
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Let Qi'^ and each stand for either / or Recalling that and that 

Qi ~ A/c(0, iV”^/) is a complex Gaussian vector, we compute 

+trQf . (5.23) 

To bound the right side of (5.23) we observe that |(QiG[®1Q2)m| ^ \g\} \ < 1, where we used 
that Bi is an eigenvector of and (5.13). (Notice that, for simplicity, here and at several 
other places we consistently use the notation ^ even when the stronger < or < relations 
would also hold, i.e., we use the concept stochastic domination even for estimating almost 
surely bounded or deterministic quantities.) 

To control the second term on the right side of (5.23), we note that a first order Neumann 
expansion of the resolvents yields 

|trg«Q«G[®l -trQ^Q^G^I = |tr G^(5, + ^^(z) - z)e,e*G« | 

l||Q«Q«G[®le,e*||2||e,e:G«||2 






^||G[®le,e:||2||e,e:G«|h 


N 


^l|(e:|G[®lre,)|'^'|(e:|G«|^e,)|'^', 


(5.24) 


where we used the boundedness of 6^, ujb{z), ||Qi*^|| and ||Q 2 *^||. Notice next the identities 

ImG'*](z) 


(|G[®l(z)n,, = 


r\ 9 I™ G^*^ (z) 

(|G«(^)n,y =-(5.25) 


V 


(1 - Sij)r] + SijlmuJBiz) ’ 

for j £ Jl, TV], with z = E + ig and |Gp = G*G. The identities in (5.25) follow directly from 
the definitions in (5.12), (4.4) and the definition of the Green function. Since |G-f | ^ ^ 
(c./., (5.5)), we obtain combining (5.24) and (5.25) with (5.3) that 

1 


|trQ«gf GW -trQ«Q«G«| ^ — . 


(5.26) 


Since iL^®^ is a Hermitian finite-rank perturbation of H, we can apply (3.3) to conclude that 


ItrQf Qf>G-trQ<*>Qf G«| 




Ng 


(5.27) 


We will now show that tr Qj'^G^®^is bounded. Using the resolvent identities and tr R^®^ = 
trR = 0, we get 

trS«G« = l-tr(^-z)G<®\ tr R<®>G<®>B<®> = z + tr {A - z)G^^\A - z), 

thus to control tr we need to bound tr {A — z)^G^®^ for = 0,1, 2. Since iV^®^ 

is a Hermitian finite-rank perturbation of H, we can apply (3.3) to conclude that 

1 


|tr {A - z)'®G« - tr (H - z)'®G| ^ 


Ng ’ 


k = 0,1,2. 


(5.28) 


Since A is diagonal with bounded matrix elements, we have tr {A — z)^G < maxj \Gjj \ -< 1, 
where the last bound comes from (5.3). This directly controls tr and then, using 

(5.26) and (5.27), we have 

|trQ<*>gf G<®>| + |tr G| + |trg<*>gf G^l ^1. 


(5.29) 
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Thus, returning to (5.23), we showed 

(5.30) 

Using the Gaussian concentration estimates in (3.1) and Wi = Bi + g^, we obtain 

'(IQfGW(z)g«|2),A" 




-< 


N 




-< 


( ImGi:'(^) ^ 

, /^ImtrGW(z)\ 

1 lVIma;B(z) , 

) +( Ng ) 


-< 


\/Nrj ’ 


(5.31) 


where we also used that is an eigenvector of that is bounded and (5.25). In the 
last step (5.13) and (5.29) were used. Combined with (5.30) we thus proved 

w*QfG^^\z)Qfw,^l. (5.32) 

For a later use we remark that, combining (5.27) and (5.31), we also proved 

w*Q^l^G^^\z)Q^^w, = (Qf>GWQ<*>),, + trQf GQ^ + • (5-33) 

In a very similar way we get, recalling that tr B = 0 and ||i3|| ^ I, that 

1 

y/N 

To deal with terms containing four or six factors of Wi in lEg. [Sj] (see (5.22)), we use the 
following rough bound. For general random variables X and Y satisfying |X|, |F| ^ 1, we 
have 


) . (5.34) 


w*B^'''>Wi = bi+ Wg.w*B^'''>Wi = bi + 0^(^ 


[XU] = lEg^ [XjIEg^ [U ] ] + lEg^ [X]Eg^ [U]] + ffig ^ [Eg^ [X]IEg^ [U ] ] . (5.35) 

In particular we have |]IEgjXy]| ^ |IIEg.X| + |IEg.U|, where we used basic properties of 
stochastic domination outlined in Subsection 3.1. 

Then, recalling the explicit form of 5^ in (5.22) and using (5.31), (5.32), (5.34) and (5.35), 
a straightforward estimate shows that |IEg.Si| ^ and claim (i) in (5.21) is thus proved. 

We next show statement (ii) of (5.21). To compute the expectation Eg. [5^], we are going 
to use the identities 


tr BG = 1 — tr (A — z)G , tr BGB = z + tr (A — z)G(A — z ), (5.36) 

that follow from {H — z)G{z) = 1 and tr A = trS = 0. Invoking assumption (5.3) we have 


N 




with k gN. Recalling further the shorthand notation mg = and from (2.12) that 

N I 

(5.37) 


I 1 


ujb 


we get from the above that 

tr G = TOffl + 0^(A^“^), 
tr ((A - z)G) = 1 + (ujb - z)ms + 0^(iV“^ ), 
tr ((A - z)^G) = wb - 2z + (wb - z)‘^m^ + 0^{N~i). 


(5.38) 
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Thus from (5.33) we obtain 

w*G'^'‘'^Wi = ma + 


Qi — UJb 


+ 0^{N-i), 


*_g<*>G-Wiy. = _ z)ma + 

= -{ujB - z)m.a + 


Qi — UJb 




Qi — UJb 




= {ujb - z) + {ujb - z)‘^ms + ■ 
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+ 0^iN-z). 


(5.39) 


tti - UJB 

Plugging (5.39) into (5.22), using the identity uja+ ujb = z — l/ma and taking the expec¬ 
tation, a straightforward computation shows that 

{uJA - bi)( 2 uJB - a^ + bi- z)ms 


1 + ^g- [“il — 


+ 0^(iV-4). 


(5.40) 


a-i — UJb 

Then from Lemma 3.4 one observes that statement (ii) of (5.21) holds. In fact, the first 
term on the right side of (5.40) is bounded away from zero uniformly on z £ 5x(?7m, !)• 

We move on to statement {in) of (5.21). Let and each stand again for either I 
or Then we note that 


= (g«GWQ«)., + (^) 


(5.41) 


as follows from the Gaussian large deviation estimates in (3.1), assumption (5.3) and the 
fact that Bi is an eigenvector of Q^i \ and G^®!. Having established (5.41), it suffices to 
recall (5.33) and (5.34) to conclude that |IEgj4'i^i]| ^ This proves claim {in) in (5.21) 

and thus completes the proof of Lemma 5.1. □ 

Corollary 5.2. Suppose that the assumptions of Theorem 2.5 are satisfied and let j > 0. 
Fix z = E + iri € 5x(?7m, 1) and assume that 

\Gif{z)-{ai-uJB{z))~'^\^N~i, \Gii{z)-{ai-ujB{z))~'^\F N~i, (5.42) 

hold for all i G NJ. Letting Q^i\stand for I or and letting Xi, stand for 
or Bi, we have the bound 

I 1 (5.43) 


max |a:*Q^*^GW(z)gWyJ ^ 1. 


i6[l.Afl 

In particular, |S'i(z)|, |T'i(2:)| ^ 1, for all i £ |1, A^]. Moreover, we have 

%jT,(z)]|^ 1 1 


max I 




max 


J5,(z)]| ^ 


■ 


(5.44) 


Proof. Using once more (3.2), we can write 


X 


*g«G(®)g«y, = =r*gfG{®}g«y^ - 


{b, +UJB-Z) a:*gf G^le, e*G{®}g« 


1 + {bi + UJB - z)gI-^ 

Hence to prove the bound in (5.43) it suffices to bound a:*gj*^G^®^g 2 ^r/j and a;*g^*^G^®^ei 
with the choices gf\ g^*^ = / or H^®^ and Xi,y^ = or Bi. To do so, we expand G^®^ 
around G^. It turns out that a;*gj*^G^®^g 2 *^yj and x*Q^^^G^^^Bi both are of the form 
4'i/(l -f Si), where is given in (5.18) and 4*^ is a polynomial of the quantities appearing 
in (5.33), (5.34) and (5.41). Then {i) and {ii) of (5.21) imply that (1 + Si)“^ ^ 1, which 
together with the bounds in (5.33) and (5.41) leads to the conclusion (5.43). 
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To prove (5.44), we follow, mutatis mutandis, the proof of (5.4) by replacing by 
Ti = or Si = For instance, for Ti the counterpart of (5.8) is 

_ gtG^^^e, - {h +u^B- z)m,^ \gI:^] gtG^'^e, 
l + {h+u:B-z)V.,XGf] 

Now, according to (5.11), (5.10) and the bound \Ti\ -< 1 (c./., (5.43)), it suffices to show 

|]DE,Jg:GWe,]H^. (5.45) 

The proof of (5.45) is nearly the same as the one of (5.10). One can also use a similar 
argument for Si by using the bound \Si\ ^ 1 from (5.43). We omit the details. □ 


6. Identification of the partial expectation Eg. [G^f] 

In this section, we estimate the partial expectation Eg^ [G\f ], which together with the 
concentration inequalities in Lemma 5.1 lead to the following lemma. Recall the definition 
of Si and Ti in (5.1). 


Proposition 6.1. Suppose that the assumptions of Theorem 2.5 are satisfied and let 7 > 0. 
Fix z = E + ig € Sxfpm., !)■ Assume that 


G^iz) - (a^-uJBiz))' 


< N- 


Gii{z) - [oi - ujb{z)) 


^ N- 


( 6 . 1 ) 


hold uniformly in i £ |l,iV|. Then, 


and 


max 

*eIi.Nl 


- (a* - Wb(^)) ^ 


1 


max 


Si{z) 


z-^%{z) ^ 1 

Qi-LUg^z) \^Nrj’ 


max \Ti{z)\ -< 
ieli.N] ' 


1 

y/Ng 


( 6 . 2 ) 


(6.3) 


In the proof of Proposition 6.1 we will need the following auxiliary lemma whose proof is 
postponed to the very end of this section. 


Lemma 6.2. Under the assumption of Proposition 6.1, the estimates 

|tr {B^^G^^{z) - BG{z)) | < ;^ , |tr G^'^{z)B^^ - BG{z)B) | < ;^ , (6.4) 

and the bounds 

|tr (i?«G6)(z))| ^ 1, |tr(il<*>GW(z)i?<*>)| ^ 1, (6.5) 

hold uniformly in i £ |1,N^]. Furthermore the estimates 

|]DEgJtr(B«GW(z))]| <^, [lEg Jtr (R«GW (z)R«)] | < ^ , (6.6) 


hold uniformly in i £ |l,iV]. 

Proof of Proposition 6.1. Fix i £ |1,A^]. By the concentration results of Lemma 5.1 and 
Corollary 5.2, it suffices to estimate Eg. [G^f (z)]. Eg. [5'i(z)] and Eg. [Ti{z)] to establish (6.2) 
and (6.3). Recall the definition of and G^^ from (4.8). We start with the identity 

{A-z)G^^Hz) =-B^^^G^^\z)+I, 


z£C+. 


(6.7) 
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Since A is diagonal, we have 

(a.-z)G«(^) = -(B«GW(^)),^ + r 


€ C+. 


( 6 . 8 ) 


Therefore, to estimate EgJG-*^(z)], it suffices to estimate instead. Re¬ 

calling the definitions in (4.6) and (4.7), we have 

(rWgW)^^ = e:(l- e^e* - e^g* - g.e* - g,gt)B^'^ 

X (/ — e^e* — Big* — g^e* — gigi)G^^^ei 

= - e *[ eig * 4- g ^ e * + g ^ g *) B ^’'^ (/ - e^e* - e ^ g * - g ^ e * - g . ig * i ) G ^'‘^ Bi . (6.9) 

Since Ci is an eigenvector of (c./., (4.3)), we have (B^'^G^®^).. = biG^f. Since moreover B 
is traceless by assumption (2.25), we have trB^*^ = trB = 0. Thus the apriori estimates 
in (6.1), the bound in (5.43), and the following concentration estimates (c./., Lemma 3.1) 


\^j9^ 


y/N ^ " y/N 


( 6 . 10 ) 


for all j € [l,iVl, imply that g*B^^>G<^'^ Bi is the only relevant term in (6.9). Thus recalling 
from definition (5.1) that Si = g*B^'^^ Bi we arrive at 


|(B«G«),,+B,| ^ 


( 6 . 11 ) 


Using integration by parts for complex Gaussian random variables, we compute Eg^ [B^] 
next. Regarding g and 'g as independent variables for computing dgf{g,g), we have 


. 1^1 


|gr 


5/(5,5) e ^d5Ad5 = cr / dgf{g,g)e ^dgAdg 


( 6 . 12 ) 


for differentiable functions / : —?► C. Using (6.12) with = 1/N for each component of 

di = {gn, ■ ■ -^gm), we have 


N 


EgJB,] = ^EgJg,,(B^>G«)fc.] = ^ E 


N 


B(B«GW)fe, 


dg. 


'ik 


(6.13) 


Using the definitions in (4.6), (4.7) and regarding gik, gik as independent variables, we have 

dW, 


dgih 


— ^kQi •) 


(6.14) 


so that 

a(B«G«) 


ki 


= bIB^^G^^ {BkB* + Bkg*)B^^ (/ - e,e* - B,g* - g,B* - g,g*)G^^^B, 


^Qik 

+ e*B«G« (/ - e^e* - - g,B* - g^g^^ [B^e* + Bkg:)G<~^^e,. (6.15) 
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Since is an eigenvector of with eigenvalue bi, we further get from (6.15) that 


- hg*G^^e,) 

+ +g*G«e,) 

- (G« +g:G«e,) {e*W'^g^+g*me^+g*mg:) 

- (G« +g*G«e,)(e*.g«efc)(e*B«G«e,) 

- (G« +g*G«e,)(g:B<*>e,)(e*B<*>G«e,) 

- (G« +g*G«eO«-§^*^‘^fc)(4S«G«9,) 

- (G« + g*G«e,) {el^^G^^g,) . (6.16) 

Plugging (6.16) into (6.13) and rearranging, we get 


EgJ^i] = Eg, [tr (B<*> G«) G«e, - b, g^G^e,) 

+ [tr (i?<*>G(')i3<*>) (gI*^ + g*G(*)e,) 


-E, 


'9i 




■tr(i?«G«) (e*il«g,+g:B«« 

'(6-G«+9*(5«)"g(' 


+e* 


+ g*B«g,)(G«+g*G«e.), 

(5«)2G«g,+g:(5«)"G«g,) 






We next claim that the last two terms on the right of (6.17) are small. Using the bounded¬ 
ness of G^f (following from the apriori estimate (6.1)), the bound (5.43), the concentration 
estimates in (6.10), and estimate (6.5) of the auxiliary Lemma 6.2, and the trivial bounds 

|^*(jj<z))2GWy^| ^ i, aji, = e, or g,, (6.18) 

y 

we see that the last two terms on the right side of (6.17) are indeed negligible, i.e., 


Eg.[^.]=Eg 


tr (B«G«)(5,-6,r,) 

+ Eg^ [tr (B<*>G(*).B«) (G« + T,)] + (^) + (^) , (6.19) 


where we also used the definitions of Ti and Si in (5.1). From assumption (6.1) and Corol¬ 
lary 5.2, we have the bounds 


max |G^'(-^|^1, max irJ^l, max ISd ^ 1. (6.20) 

J6I1WI " teliWI jelbJVl 


We hence obtain from (6.19), (6.5), and the concentration estimates in (6.6), (5.4) that 
Eg^ [5,] = tr (B<*>G(*)) (Eg^ [5,] - h Eg, [T,]) 

+ tr(B«G«5«) (G« +EgJr,]) +0^(^) . (6.21) 
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Repeating the above computations for Eg. [g* Cj] = Eg. \Ti], we similarly obtain 
Eg^ [T,] = tr G(*) (Eg^ [5,] - h Eg, [tJ ) 

+ tr (R<*>G(*)) (G« + Eg^ [T,]) + (^) . (6.22) 

Now, using the bounds in (6.20), the estimates (6.4) and |trG^*) — trG| ^ (following 

from (3.3)), we obtain from (6.21) and (6.22) the equations 

Eg^ - tr (RGR) (G« + Eg^ [T,] ) = tr {BG) (Eg, [5.] - 6.Eg, [T,] ) + (^) , 

(6.23) 

and 

Eg, [T,] - tr [BG) (G« + Eg, [T,] ) = tr (G) (Eg, [5,] - 6.Eg, [T,] ) + (^) . (6.24) 

We first approximately solve (6.24) for Eg. [T^] to show, under the assumptions of Propo¬ 
sition 6.1, that jEg.Til ^ N~i. To see this, we recall (6.8) and (6.11) which together with 
assumption (6.1) imply that 

Si = {tti — 2)G-j^ -1-1-0^ (“7^) ~-h (A^ <1) . (6.25) 

N ' di ~ dJB 

By the concentration estimate (5.44), we also have 


Eg,[5,] =- 


z-lub 


+ 0^ {N-i). 


(6.26) 


di — U!b 

In addition, by the identity BG = I — {A — z)G, assumption (6.1) and equality (5.37), we 
have, using the shorthand notation 

trG = TTiffl -I- (A^“^) , tr (^BG) = {z — ujB)m^ + . (6.27) 

Substituting (6.26) and assumption (6.1) into (6.24), and using \Ti\, jS'^j ^ 1, we obtain 


|(l-tr(RG) -f 6,trG)Eg, [T,] | ^ N-i . 

Using (6.27) and the second equation of (2.12), we have 

I (l - tr {BG) + bitr G) | = 11 -I- {ujb - z + bi)m^ | -|- (-^"^) 

= \{-ujA + bi)ms\ + {N~^) . 

Since |(—+ bi)m^\ > 1 by (3.15), we have from (6.28) that Eg, [T^] ^ N~i. 
from (5.44), [Til ^ N~i. Then solving (6.23) and (6.24) for Eg, [^i], we obtain 


(6.28) 


(6.29) 

Hence 


lEg, [S^] = Gf) + 


trG 
-f 1 


tr {BG) - (tr (HG))" 




trG 


-tr {BGB) (G«+Eg,[r,]) 


(6.30) 


Averaging over the index i and reorganizing, we get 


tr(BG) - (tr {BG)y 


trG 


+ tr {BGB) 


+E..M) +0^(^) 


Nrj ^ 


^E"i(Gl:^+EgjT,]) 


(6.31) 
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Now, recalling the concentration of Si in (5.44) and estimate (6.11), we have 

1 




y/Nr] 


(6.32) 


Note that under assumption (6.1), we can use Corollary 5.2 to get (5.43), which together 
with (6.1) implies that the assumptions in Lemma 4.1 in the case oii = j = k are satisfied. 
Then, by (4.10) with i = j = k and (6.8), we get 

1 




o. 


for all i G |l,iV]. Using (6.32) and (6.33) we obtain 




N 


^G«-trG 




1 


i//V?7 ’ 


N 


[5.] +tr(BG) 


1 




y/W^ ' 


(6.33) 


(6.34) 


'' ' i=l 

Substituting (6.34) and assumption (6.1) into the right side of (6.31), and using |trG| > 1 
(following from (6.27)) and \Ti\ -< N~i, we obtain 

> 2 


tr {BG) - (tr (BG))" 


trG 


Now, plugging (6.35) back into (6.30) gives 


+ tr {BGB) 




1 


Eg. 


trG 


//V^ ■ 


which together with (6.8) and (6.32) implies that 


(6.35) 


(6.36) 


(6.37) 


in light of the definition of uj%{z) in (2.28). By assumption (6.1) we see that u!g{z) = 
u!Biz) + 0^{N~i). Hence by (3.15), we also have lmujg{z) > c for some positive constant c. 
Therefore, we get (6.2) from (6.37). 

Then (6.36) and (6.2), together with the definition of ujg{z) in (2.28) and the concentra¬ 
tion of Si in (5.2), imply the estimate of Si in (6.3). 

Substituting (6.36) into (6.24), we strengthen (6.28) to 

I (1 - tr (BG) + 6,tr G)Eg^ [T,] | ^ . (6.38) 


1 


Using (6.29) again, we obtain from (6.38) that 

lEg, [Tr] I ^ ^ , 

which together with the concentration inequality in (5.44) implies (6.3). Therefore, we 
complete the proof of Lemma 6.1. □ 

We conclude this section with the proof of Lemma 6.2. 

Proof of Lemma 6.2. We start by invoking the finite-rank perturbation formula (3.3) to get 

|trQf GWQ^*^ - trQf^GQ^*^! < , 


=/orB«. 
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Hence, it suffices to verify (6.4) and (6.5) with replaced by G. Recalling from Section 4 
that Ri = I — riV* and using the fact that Ri is a Householder reflection (in fact \\ri \\2 = 2 
by construction), we have = RiBRi. Then we write 

tr (B<®>G) = tr {R,BR,G) = tr (BG) + di , (6.39) 

with 

d, := -^r*BGri - ^r*GBr, + j^{r*Br,){r,Gr,). 

Using that ||G|| < l/r], we immediately get the deterministic bound \di\ < G/Nr], for some 
numerical constant G. Together with (6.39) this implies the first estimate in (6.4). The 
second estimate in (6.4) is obtained in the similar way. 

The bounds in (6.5) follow by combining the sharp formulas for tr (BG) and tr(BGB) 
from (6.27), (6.35) with the estimates in (6.4). 

To prove ( 6 . 6 ), we set = B^®^ or (B^®^)^ and note that 

jlEgJtr (Q«G«)]| = [iDEgJtr (Q«(G« - G«))]| < , 

where we used that and G^®^ are independent, and once more (3.3). □ 

7. Proof of Theorem 2.5: Inequalities (2.19) and (2.21) 

In this section, we prove the estimates (2.19) and (2.21) of Theorem 2.5 via a continuity 
argument. We also prove Theorem 2.6. 

First, let us recall the matrix R and its Green function Q defined in (2.26) and (2.27), 
these are the natural counterparts of H and G with the roles of A and B interchanged. We 
can apply a similar partial randomness decomposition to the unitary U* in R as we did 
for C/ in B in Section 4. This means that, for any i £ |l,iV|, there exists an independent 
pair {vi,W), uniformly distributed on and U{N — 1), respectively, such that with 

ri := \/2(ei + e~^^'Vi)/\\ei + e“'®®9i||2, we have the decomposition U* = —e'®®77.iG^®\ 
where 6i is the argument of the ith coordinate of Vi] where Ri := (I — rir*) and is 
the unitary matrix with as its ith column and G® as its (i,i)-matrix minor. Analogously 
to Qi defined in (4.5), we define a Gaussian vector g, = (g^i,... .gm) ^ A/c(0, B“^/), to 
approximate Setting Wi := and Wi ■= I — WiW*, we dehne 

B^®) := B +W,G<®>A(G<®>)*>Vi, 

for all i G |1, N}. Calligraphic letters are used to distinguish the decompositions of R from 
the decompositions of H. 

Next, we introduce the z-dependent random variable 

Ad(z) ;= max |G-)^(z) - (oi - 015 ( 2 )) ^| + max |Gij(z) - (a* - 015 ( 2 )) 

+ max \gll\z) - {b, - uja{z)) ^|+ max \g,i{z) - {b, - uja{z)) Y (7-1) 

iell.N] i 6 [l.Af] 

Moreover, for any 6 G [0,1] and 2 G Si{r]m, 1), we define the following event 

0d(B<5):={Ad(2)<^}. (7.2) 

The subscript d refers to “diagonal” matrix elements. With the above notation, we have 
the following lemma. 
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Lemma 7.1. Suppose that the assumptions of Theorem 2.5 are satisfied and fix 7 > 0. 
For any e with 0 < £ < ^ and for any D > 0 there exists a positive integer N 2 {D,e) 
such that the following holds: For any fixed z = E + \r] £ Si^rj^n, 1) there exists an event 
ftd{z) = fldiz,D,£) with 

¥{nd{z))>l-N-^, V 7 V> 7 V 2 (Ae), 

such that if the estimate 

F{ed{z,N-i)) + (7.3) 

holds for all D > 0 and N > iVi(Zl, 7 ,e), for some threshold Ni{D.,j,£), then we also have 

ed{z,N 4)niid(z) c (7.4) 

for all N > N 3 {D,j,£) := max{iVi(Zl, 7 , e), A^ 2 (L^, e)}- 


Proof. In this proof we fix 2 ; € Sx{r]m, !)• 
that assumption (7.3) implies 

- K - ^ > 

and 

- ^a{z))~^\ -< , 


By the definition of ^ in Definition 1.1, we see 
\Gu{z)-{a,-u;'fi{z))-^\^N-i , (7.5) 

\g^^{z)-{h-UJ^A^))-^^N-i. (7.6) 


Hence, we can use Corollary 5.2 to get (5.43). Together with the boundedness of 
and Gii {c.f, (7.5) and (3.15)) this implies that the assumptions in (4.9) of Lemma 4.1 are 
satisfied when i = j = k. Thus (4.10) holds when i = j = k. Hence, invoking, (7.5) and 
Proposition 6.1, we get 

|G«(z) - (a. - ^ , \Gu(z) - (a. - -< ^ ■ (7.7) 

Switching the roles of A and B as well as U and U*, and further using (2.29), we also get 

\gjp(z) - ih - u:\(z))-^\ F ^ , \gu(z) - ih - uj'Xiz))-^ F ^ , (7.8) 


under (7.6). 

Now, we state the conclusions (7.7) and (7.8) in a more explicit quantitative form assum¬ 
ing (7.3) which is a quantitative form of (7.5)-(7.6). Namely, we show that the inequalities 


T\T — 

G\f(z)-(a.-u;Uz))-^\<^, 

|t;«(z)-(6,-.;i(z))-i|<^. 


Guiz) - (a, - uj'fi(z)) 

\gu{z)-{bi-uj'x(z))~^\<-^=. (7.9) 


hold on the event Qd(z, N~i) fl Dd(2:), when N > N 3 (D,"f,£). Here Dd(z) is an event 
determined as the intersection of the “typical” events in all the concentration estimates in 
Sections 4-6. 

To see this more precisely, we go back to the proofs in these sections. The concentration 
estimates always involved quantities of the form W,g\g*Qx] with x = gi,e.i and some 
explicit matrix Q that is independent of g^ but often z-dependent. The total number of 
such estimates was linear in TV. Thus, according to Lemma 3.1, for any (small) £ > 0 and 
(large) D > 0, there exists an event gid(z,D,£) with 


P(Dd(z,D,£)) > 1-iV-^ 


(7.10) 
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such that all estimates of the form 

|IEgJg:Qe,]|<^||Qe,|| 2 , \m,^[g:Qg,]\ < ^\\Qh (7.11) 

in Sections 4-6 hold on D, e) for all N > N 2 {D, e). In addition, the threshold N 2 {D, e) 
is independent of the spectral parameter z. 

We now follow the proofs in Sections 4-6 to the letter but we use (7.10), (7.11) and (7.3) 
instead of the ^ relation. Instead of (7.7) and (7.8), we find that the analogous but more 
quantitative bounds (7.9) hold on the intersection of the events 0d(z, N~i) and fld(.z, D, e). 
It remains to show that on the event 0d(^, N~i) fl nd(.z), 

CNi CNi 

- u;a{z)\ < , l^l(z) - u^b{z)\ < (7.12) 

hold when N > iV 3 (_D, 7 ,e). 

To this end, we use the stability of the system ^z) = 0 as formulated in 
Lemma 3.3. By the dehnition of the approximate subordination functions and ^^(z) 

in (2.28), by the identity (2.30) and by taking the average over the index i in the estimates 
in (7.9), we get the system of equations 


niHiz) = mA{u}B{z)) + rA{z) , 
mi/(z) = TOB(a;^(z)) -|- rsiz) , 


a;^(z) -hWB(z) = z- 


mniz) 


where the error terms ta and tb satisfy 


\rB{z)\ < 


CNi 

7^’ 


(7.13) 


on the event Qd{z,N~i) fl nd(.z) when N > Using the definition of 0d(z,(f) 

in (7.2), (7.9) and the fact that z G 5x(?7m, 1), so oja{z) and u}b{z) are well separated from 
the real axis, we have 


\uj^Az)-^A{z)\<CN-i , \u;Uz)-uJB{z)\<CN-i . (7.14) 


on the event 0d(z, N 4 ) nUd(z) when N > N 3 {D, 7 , e). Hence, plugging the third equation 
of (7.13) into the first two and using (3.15) together with (7.14), we get 

^ma.mb(wa(^:),Wb(z),z) =Hz), 

where r(z) = (rA(z), rB(z)) with 

CNi CNi 

P.15) 

on the event &d{z,N~^) fl Hd( 2 :) when N > A^ 3 (_D, 7 ,e). Therefore, by Lemma 3.3, we 
get (7.12). Hence, we completed the proof of Lemma 7.1. □ 


Given Lemma 7.1, we next prove Theorem 2.5 via a continuity argument similarly to [19]. 


Proof of (2.19) of Theorem 2.5. From Theorem 1.2 of [26], we see that for 77 = 1, we have 

max \Cii{z) - {ai - u}Biz))~^\ ^ N~i, max \gi^{z) - {bi - u!Aiz))~^\ P N~i 
ieli.TVl I ' ^ ^ I iGlpJV] I ' ^ ^ wy I 


(7.16) 
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if 0 < 7 < 1/7 (say). In addition, owing to the estimate ||G|| < l/ry, assumption (4.9) 
obviously holds for 77 = 1. Hence, by Lemma 4.1 in the case oi i = j = k and its analogue 
for we have 


max IG(z) - (a^ 
zGll.W]' 


usiz)) M ^ IV 2 ^ 


max \Q\f{z) - {hi - uja{z)) M ^ IV 2 . 

(7.17) 


Hence, for any E €l and D > 0, 

F{Qd{E + i,N-^))>l-N-^, (7.18) 

holds for all N > No{D,'-f) with some IVo(Il, 7 ) > 0. In the sequel we will apply Lemma 7.1 
with the choice 


Ni{D,'y,e) := max{No{D,j),N 2 {D,e)} ; 

in particular we have N 3 {D, 7 , e) = Ni{D, 7 , e). 

Next, we define the lattice 

1) := 5 i( 77 „„ 1) n IV^jZ x iZ} . (7.19) 

Thanks to the Lipschitz continuity of the Green function, i.e., ||G( 2 ;) — G( 2 :')|| < IV^|z — z'\ 
for any z,z' G Sz(r)m,l), it suffices to show (2.19) on the lattice 5 i( 77 ni,l). We now fix 
E G Tn IV“®Z and decrease 77 from 77 = 1 down to steps of size IV“®. Recall the 

events 0d(2, <5) and Hd( 2 ;) in Lemma 7.1, and choose the same e < ^ in fid(z,E,e) for all 
z. For simplicity, we omit the real part E from the notation and rewrite 

0d(77,S) := Od(E + i77,<5), Hd(77) := fid(^ + iy)- 

Our aim is to show that for any 77 G [ 77 m, 1], 

0d(77,IV-T) nHd(r7-IV-®) C 0d(77-IV-5,IV-T). (7.20) 

To see (7.20), we first notice that by the Lipschitz continuity of the Green function and of 
the subordination functions oja(z) and ujb(z) (see (3.9)), we have 

0 d( 77 ,IV-^) C 0 d( 77 -IV-^Af-^ +GIV-3) c 0 d( 77 -IV-^Af-^), (7.21) 

where the last step is obtained by choosing 7 > 0 sufficiently small. Now, we start from 
(7.18). By (7.21), we get 

P(0d(l - IV-^IV-^)) > P(0d(l,IV-^)) > 1 - iv-^. 

Hence, we can use Lemma 7.1 to get 

e,(i - JV-=. «-}) n na(i - c e, (i - iv-=, 7 ^ 7 ^=^) 

C ed(l - , (7.22) 

which together with (7.21) implies (7.20) with 77 = 1. Now, replacing 1 by 1 — IV“®, we get 
from (7.22), (7.18) and the fact P(Hd(l — IV“®)) > 1 — for IV > N 2 iD,e) that 

P(0d(l - IV-^IV-^)) > 1 - 2^--° (7.23) 

holds for all N > N^{D,^,e). Now, using (7.23) instead of (7.18), we get (7.20) for 77 = 
1 — IV“^. Iterating this argument, we obtain for any 77 G [rjm, 1] H N~^'L that 

0 d(l,IV-T) n Hd(l - N-^) n ... n Hd(r7) c 0 d( 77 , . 
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Hence, we have 

P(0d(77,A^”^)) > l-iV-^(l + A^^(l-r?)) 
for all N > A^ 3 (D, 7 ,e), which further implies 

P(0d(?7-iV■^A^■^)) > 1- N-^{l + N^l-r])) 
for all N > A^ 3 (£), 7 ,e), by using (7.21). Then, using Lemma 7.1 again, we obtain 

P(0d(r;, ^)) > 1 - iV -^(2 + N^l - rj)) (7.24) 

uniformly for all rj G [?7m, 1] H N~^h, when N > A^ 3 (Il, 7 ,e). Finally, by continuity, we can 
extend the bounds from z in the discrete lattice to the entire domain 5x(7m, !)■ We then get 


max 

ielpivi 


Gu{z) - 


1 


a* - ujb(z) 




1 


max 

ielb^l 


G^{z) - 


1 


ai-UJB(z) ^/Wq 


-< 


1 


(7.25) 


v^’ 

uniformly on 5x(7m,l), where we used the definitions of 0 d( 2 :,< 5 ) in (7.2) and of ^ in 
Definition 1.1. This concludes the proof of (2.19). □ 

Having established (2.19), (2.21) of Theorem 2.5 and Theorem 2.5 are direct consequences. 

Proof of (2.21) of Theorem 2.5. It suffices to note that (2.21) is a direct consequence of (2.19) 
and the facts mniz) = Gu{z) and mASB{z) = iV“^ “ wb)“P □ 

Proof of Theorem 2.6. Using the spectral decomposition of the Green function G, we have 

max lTa.Gjj{z)= max 'Y' ———2 = 2 € C+. (7.26) 

icn Aril ^ ' ocn i\ni ^ X.- — K 2 1 „2 _ Tr \2 1 „2 ’ ^ / 


jell.Afl 


|Ai - Fip + 772 ^ |y. _ £;|2 ^2 ’ 


Fix a small 7 > 0. For any Ai G X, we set U = on the right side of (7.26) and use (2.19) 
to bound the left side of it with z = \i + \q^ q = 7V“^+'''. Then we obtain 

II«*IIL ^ ?7 = iV^+^. 

Since 7 > 0 is arbitrarily small, we get (2.24). This completes the proof of Theorem 2.6. □ 

8. Proof of Theorem 2.5: Inequalities (2.20) 

In this section, we prove (2.20) of Theorem 2.5. Note that, from (7.25) in the proof 
of (2.19) in Theorem 2.5, we know that the following estimates hold uniformly on Sxfq-m, 1), 


G^h^) - {ai-ujB{.z)Y 


-< 


1 


^/Nq ’ 


Gii{z) - (oi - iJB{z)) 


-< 


1 


^/Nq ■ 


( 8 . 1 ) 


Taking (8.1) as an input, we follow the discussion in Sections 5-7 to prove the estimate (2.20) 
with the following modifications. We introduce the quantities 

X,, (z):=g:GW(z)e,, 5,,,(z) := g*B«GW(z)e, , z G C+. (8.2) 

that generalize Ti[z) and Si[z) defined in (5.1). In particular, Ti[z) = Ti^i{z) and Si{z) = 
Si^iz), but we henceforth implicitly assume that i ^ j. (We use a comma in the subscripts 
of Sij since they are not the entries of some matrix.) We often abbreviate Tij = Tij{z) 
and Sij = Sij{z). 

We first establish the concentration estimates for Gfj (see Lemma 8.1), and Tij and Sij; 
see Lemma 8.2. In Proposition 8.3 we then derive self-consistent equations for Eg^ [I'ij] s-nd 

Eg. that will show, together with concentration estimates, that |G^*^|, \Tij\, ^ 

kO 


/Nr] 


, provided that |G)!- | ^ 1. We then close the argument via continuity. 
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We start with the analogue of Lemma 5.1 for the off-diagonal entries of 


Lemma 8.1. Suppose that the assumptions of Theorem 2.5 are satisfied and let 7 > 0. Fix 
z = E + IT] G Sx{r]m, 1) and assume that 


for all i,j G |l,iV], i ^ j. Then 




(8.3) 

(8.4) 


\/Nri ’ 

for all i,j G [l,iV|, i ^ j. 

Proof of Lemma 8.1. Recall and defined in (5.12) and (5.6), as well as their Green 
functions G^ and G^®^. In the proof of Lemma 5.1 we derived the identity 


= 


Gg^(^) - {h+^B{z) - z)m^}G\l^{z)] Gt’iz) 


w. 




(8.5) 


l + {bi+ ujb(.z) - zfEg^ [GL (z)] 

see (5.8). With (8.1), we see that assumption (5.3) is satisfied. Hence, we can use all 
the conclusions in the proof of Lemma 5.1. Therefore, according to (5.11), (5.10) and 
assumption (8.3), it suffices to show the concentration estimate |IEgJG|^*^(z)]| ^ 
this end we expand G^®^ around G^®!. Recall from (5.17) that 


To 


^{*1 _^[»] 


- 


1-f: 


where :z.i is defined in (5.18) and 'i'ij is defined in (5.19). Recalling statements (i) and (ii) 
in (5.21), it suffices to establish that lIEgJ'I'ij]| ^ Note that 'I',;.,- contains the terms 

listed in (5.33), (5.34) and (5.41), as well as the terms 


e*_B«G[®lej, 


1*G^^t 


S'* 


_B<®)g[®1< 


( 8 . 6 ) 


Since is an eigenvector of B^®^ and of BW, we have e*B^®^G[®lej = SijbiG^^. Moreover, 
using Lemma 3.1 with Q^®^ = / or B<®> , we have 


/ImG 


(g<®>G[®le,| ^ _||g«G[®le,|U ^ 


(8.7) 


To control G^*j, we recall from (5.12) that the matrix ijl®! is block-diagonal and we thus 
have, for j ^ i 

Gg = (H® + WB^iur - zlN-i)il , i := jHj < *) + (j - l)l(j > ^), 

where H® and B® the are (i, i)-matrix minors of A and B respectively (obtained by removing 
the ith column and ith row) and B® G U{N — 1) is the (i,i)-matrix minor of t/^®^ which is 
Haar distributed as seen at the beginning of Section 4. Note that the matrix H® -|- [/®B®([/®)* 
satisfies the assumptions of Theorem 2.5. We thus have the estimate 


max I G^f] (z) I ^ 1. 


Z e 5i(r7m, 1). 


( 8 . 8 ) 


Plugging this bound into (8.7) we get |g*g^®^G[®lej| ^ 


/TGrj' 


The remaining part of the 


proof is nearly the same as the one of Lemma 5.1. We omit the details. 


□ 


We have the following analogue of Corollary 5.2. 
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Lemma 8.2. Suppose that the assumptions of Theorem 2.5 are satisfied and let ■j > 0. For 
all z = E + ir] € 5x(?7m, 1), we have the bounds 


max \Tij{z) \ -< 1. 


and the concentration estimates 

max llEg. [Ti^j(z)]I ^ 




max |S'jj( 2 ;)| ^ 1 


max llEg. [S'i^j(2:)]| ^ 
i^j I » ’ 


\/Nr] 


(8.9) 


( 8 . 10 ) 


Proof. With the estimates in (8.1) and (8.8), the proof is analogous to that of Corollary 5.2. 
Here we get the conclusions for aX\ z = E -\-ip £ iSx(?7m, 1) at once, since we use the uniform 
estimate (8.1) instead of assumption (5.42) for one fixed z. We omit the details. □ 

Finally, we have the following counterpart to Proposition 6.1. 

Proposition 8.3. Fix z = E + ip £ Sx{prm !)• Under the assumptions of Lemma 8.1, we 
have 

max |G'-*^( 2 :)| ^ 


and 








( 8 . 11 ) 


( 8 . 12 ) 


Proof. The proof is similar to that of Proposition 6.1. Having established the concentration 
inequalities in (8.4), it suffices to estimate Eg. [cfj] to prove (8.11). We then start with 

(a, - z)G[f = -(.§(*)GW)y + . (8.13) 

Choosing henceforth i ^ j, mimicking the reasoning from (6.9) to (6.11) and using (8.9), we 
arrive at 


(hWg(*)),, = (^) = -5.., + 

Then, instead of (6.17), we obtain 

Eg^ [,5,,,] = Eg^ ftr - b, T,,,) 


(8.14) 


+ Eq 


- E„ 


tr (H«G«H«)(Gg^ +r,.j) 




tr(.g«GW) +T,,,) 

(6?G« +g*(H«)"G«e, + e*(H«)2GWg, + g*(H«)"GWg,) 
X iG^+T,,l) 


(8.15) 


where we directly used the definitions in (8.2). Then, similarly to (6.23), using the concen¬ 
tration estimates in Lemma 8.1 and in Lemma 8.2, as well as the Gaussian concentration 
estimates in (6.10), the bound (6.18) and Lemma 6.2 for tracial quantities, we obtain 

Eg^ [5,.,] - tr {BGB) (G« + Eg^ [T,,,]) = tr {BG) (Eg, [5,., - 6.Eg^ [T,,,]]) + (^) . 

(8.16) 
















32 


Analogously, we also have 

Eg. [%,j\ - tr {BG) (G« + [T„-] ) = tr (G) (E^^ [S,,,] - [T,.,] ) + (^) . 

(8.17) 


Solving Eg. [Sij] from (8.16) and (8.17), we have 

■tr(i?G) - (tr {BG)y 




trG 


trG 


tr {BGB) 


(GW+EgjT,,,]) 


O- 


( \/Nr^ 


Using (6.35), the assumption |G-*)| ^ 1 and the bound \Tij\ ^ 1 of (8.9), we have 




which together with (8.13), (8.14), the concentration estimate (8.10) implies that 




(8.18) 


(8.19) 


This proves the estimate in (8.11). 

Next, we bound Si^j. Starting from (8.18) we directly get the second estimates in (8.12) 
from the Green function bound (8.11) and the concentration estimate (8.10). 

It remains to estimate Tij. Plugging the bound on Gy in (8.11) and the bound on Sij 
in (8.12) into the equation (8.17), we obtain 

I (1 - tr {BG) + htr G)Eg^ [T,,,] | ^ . (8.20) 

Invoking the estimate (6.29) we get |Eg^ [Ti^-] | ^ Then the first estimate in (8.12) 

follows from the concentration estimate for Tij in (8.10). This completes the proof. □ 


Having established Lemma 8.1 and Proposition 8.3, we next prove (2.20) of Theorem 2.5 
via a continuity argument similar to the proof of (2.19). 


Proof of (2.20) of Theorem 2.5. Fixing any z € Sx{r]m, 1) and using Proposition 8.3, under 
the assumption 

max|G)]^(z)H 1, (8.21) 

we have 

Then, by (8.22) and (8.1), we can use Lemma 4.1 to get 

|Gy(z)| ^ (8.23) 

Hence, in principle, it suffices to conduct a continuity argument from 77 = 1 to 77 = r]^ (similar 
to the proof of (2.19) of Theorem 2.5) to show that the bound (8.21) holds uniformly for 
z € 5’i(77m, 1). However, in order to show that (8.23) also holds uniformly for z G Sx{r]^, 1) 
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quantitatively, we monitor Gij in the continuity argument as well. To this end, we introduce 
the ^-dependent random variable 

Ao = Ao(z) := max \G^^hz)\ + max |Gy (z)|, 
and, for any d G [ 0 , 1 ] and z G Sx{r]m, 1 ), we define the event 


Qo{z,6) := {Ao(z) < (5}; 


c.f., (7.1) and (7.2). The subscript o refers to “off-diagonal”. 

We will mimic the proof of (2.19). Analogously, using Lemma 4.1 and Proposition 8.3, one 
shows that there exists an event no(-z) = ^o{z, D, e) such that the conclusions in Lemma 7.1 
still hold when we replace 0d(z, d) by © 0 ( 2 :, (5), fld( 2 :) by no(- 2 ) and N~4 by 1. We also set 
(5 = 1 in this proof. This is a quantitative description of the derivation of the first bound 
in (8.22) and (8.23) from (8.21). The main difference is that here 110 ( 2 :) is the event defined 
as the intersection of the “typical” events in all the concentration estimates in Sections 4-6, 
in the proofs of Lemma 8.1 and Proposition 8.3, and the event on which the following bounds 
hold 


Guiz)-{a,-u}Biz)) ^ 
Gii(z) - (a^ - ujb(z)) ^ 


< 


< 


\/Nrj 

7 ^’ 


max|G'^(z)| < . 


(8.24) 


Note that, by (8.1) and ( 8 . 8 ), we know that (8.24) holds with high probability uniformly on 

With the analogue of Lemma 7.1 for €>o(z,d = 1) and llo( 2 :), we conduct a continuity 
argument similar to the one in the proof of (2.19). Again, by Lipschitz continuity of the 
Green function it suffices to show estimate (2.20) on the lattice Sxirjim 1) defined in (7.19). 
We fix G G ZnN ~^ Z , write z = E + i-rj and decrease 77 from 77 = 1 down to in steps of 

size N~^. The initial estimate for 77 = 1, i.e., Ao(G -|- i) < 1 follows directly from the trivial 
fact ||G^®)(z)||, ||G(z)|| < 1 / 77 . Then one can show step by step that for any 77 G [ 77 m, 1], say, 

© 0 ( 77 , 1) n iloiv - N-^) C © 0(77 - 1 ) , (8-25) 


which is the analogue of (7.20). The remaining proof is nearly the same as the counterpart 
in the proof of of (2.19). We thus omit the details. □ 


Appendix A. Orthogonal case 

In this appendix, we show that Theorem 2.5 and Theorem 2.6 also hold in the orthogonal 
setup where U is Haar distributed on the orthogonal group 0{N). From the proof of 
Theorem 2.6, we see that it is implied by Theorem 2.5. Hence, it suffices to discuss the 
latter. We outline the necessary changes in the discussion of Sections 4-8 to adapt our proof 
to the orthogonal case. We mainly show the modification for the proof of (2.19) in detail, 
and ( 2 . 20 ) will be discussed briefly at the end. 

First, we modify some notation. We start with the decomposition for the Haar measure 
on the orthogonal group analogous to (2.35). For all i G |1, A^l, according to [28], there exist 
a random vector Vi = {vn,... ,ViN)^ uniformly distributed on the real unit {N — l)-sphere 
'■= {x G : x*x = 1}, and a Haar distributed orthogonal matrix U’' G 0{N — 1), 
which is independent of Vi, such that one has the decomposition 

U =- sgn{vii){I - rir*)U^"^ := - sgn{vii)Ri , 








34 


where 


h := 


■ sg'n{vu)vi 


—- , , „ , R,-.= I-r,r*, (A.l) 

l|ei +sgn(?;ji)t!i||2 

and is the orthogonal matrix with as its ith column and W as its (i, i)-matrix minor. 
Moreover, there is a real Gaussian vector g^ ^ A/r( 0, such that 


\\9i\\2 


Similarly to (4.5), we define 

gik ■■= sgn{vii)gik , k^i, 

and introduce an 7V(0,iV“^) variable ga^ which is independent of the orthogonal matrix U 
and of g^. Let g^ := {gn,... ^gm) and note that g^ ^ A/r( 0, Then we set Wi := 

ei + gfj and Wi := I — WiW* as before. With these modifications, we follow the proofs in 
Sections 4-7 verbatim. The only difference is the derivation of (6.19). Instead of (6.12), we 
use the following integration by parts formula for real Gaussian random variables 


/ 5 /( 5 ) e ^dg = a‘^ / ,/'(g)e ^dg, 

JR Jr 

for differentiable functions / : E —>■ R. Correspondingly, instead of (6.14), we have 

dW,, 4. 4: 4. 

^kQi Qi^k ' 


(A.2) 


d^gik 


Hence, we get 


kj 


dgtk 

= {e,el + eue* + etg* + g,el)B^^^ (/ - e,e* - e,g* - g,e* - g,g*)G^^h 

+ (/ - - e,g* - g,e* - g,g*)&'^ (e,e* + eke* + e^g* + 9,e*)G(*)e 


instead of (6.15). Substitution into the identity 


N 


1 


N 






dgik 


yields 


EgJg*i?«GWe,] = (r.h.s. of (8.15)) + -E^^ [(G«(i?«)2G«)^.^ 


+ 

+ 


g*GW(B«)2GWg^ 


hje*G^'^B^^G^^g, 


- —E„ 

N 

——Eq 

N 


(G(i)^<i)Gb)5(b)^ 


(G« + g*GWe,)(e*(H«)2G«g, + 

{(^ +g*G^ei){fiG^ +g*{B^fG^, 


(A.3) 


where we introduced bi := WiB^'^'^Wi. Note that the last two terms were discussed in the 
unitary setup, and they were shown to be negligible. Therefore, to get (8.16) also in the 
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orthogonal case, we rely on the following lemma to discard the supplementary small terms 
in (A.3). At first, let us discuss the case of i = j, which suffices for the proof of (2.19). 


Lemma A.l. Under the assumption of Proposition 6.1, we have the following bounds 




|e*G«(z)B«G«(z)g,| ^ i, 
for all i € |1,..., iV|. 

Proof. We drop z from the notation. For the first two terms, we have 


(GW(z)(i?«)2GW(^))^J^i, 
e*GW(z)(5«)2GW(z)g,| ^ i 


(A.4) 




Im G, 


(d 




1 


(A.5) 


V V 

for fc = 1, 2, where in the last step we used assumption (6.1). For the third and fourth terms, 
we have, for fc = 1, 2, 

|e*G«(5«)'^'GWg,,| < ||S«f ||G«e,||2|lGWgJ|2 < (((GW)*G«)^J^(g*(G«rG«gO^ 

= i(lmG«)^(lmg*G«g,)^ 

where in the last step we used assumption (6.1) and estimate (5.43). For the fifth term we 
note that = bi{G^^'lB^^'>G<^^'>)u and the bound follows from (A.5). For 

the last term, we have 


1 


where we used assumption (6.1) and estimate (5.43). This completes the proof. 


2*GW||2||G«i?«g,|l2^- 


(A.6) 

□ 


All the other arguments in Sections 4-7 work for the orthogonal case as well without 
modifications. This proves (2.19) of Theorem 2.5 for the Haar orthogonal case. 

For (2.20), analogously to (A.4), we shall estimate the second to the seventh terms on 
the right side of (A.3), under the assumption of Proposition 8.3. To bound these terms, we 
can pursue the discussion from (A.5)-(A.6) with ImG^*^ replaced by ImG^*^ in the bounding 
procedure. Hence, it suffices to show for all i that |Gj*^ 1^1- To see this, we use the fact 
\Gjj \ -< 1 from Theorem 2.5 and \Gjj — Gj*^| ^ "7^ from Lemma 4.1 with k = j. Note that 
the assumptions of Lemma 4.1 are guaranteed by Theorem 2.5, (8.9) and assumption (8.3). 
Hence, (2.20) also holds in the orthogonal case. 

Appendix B. Two point mass case 

In this section, we present our result when both, pa and pp, are convex combinations of 
two point masses. Without loss of generality (up to shifting and scaling), we may assume 
that Pa and pp are of the following form, 

Ma = + (1 - C)<5o , M/3 = + (1 - C)<5o , 

with real parameters A C and 9 satisfying 


0AO, 


e,ce (o, 


e<c, 


{^,^,0 7^ ( 


- 1 , 


1 1 

2 ’ 2 


(B.l) 
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Here we excluded the case = (— 1 , since it is equivalent to (0,^, C) = (1) 51 \) 

under a shifting, where the latter is a special case oi iJ,a = In Section 7 of [1], we 
explained why the setting of (B.l) is special, and we thus excluded it from Theorem 2.5. 
Following [24], we argued in Lemma 7.1 of [1] that in the setting of (B.l) we have 

=(^ 1 :^ 2 ) U (4,4), (B.2) 

in case ^ while we have 

= (4,4), (B.3) 

in case /Tq, = /i/ 3 , where 

h := min { ^ (l + ^ “ V (1 - 9Y + 40r+^ , i (^1 + 0 - yj{I - 9^ + 40r_^ | , 

4 := max { ^ (l + ^ “ \/(l - Oy + 40r+^ ’ ^ ^ ~ V (1 - 0 )^ + 40r_^ | , 

4 := min { ^ (l + ^ + \/(l - GY + 40r+^ ’ ^ ^ ■\/(l “ + 40r_^ | , 

4 := max { ^ (l + ^ + V(1 - + 40r+^ ’ ^ ^ a/( 1 “ + 40r_^ | , 

with r± := ^ + C - 2 ^C ± 4^4^41-4(1-0- 

In Remark 7.2 of [1] we argued that, in the case = M/ 3 , tli® point E = 1 G 
is special in the sense that m,^^[jj^^(l + iO) is unstable under small perturbations. We thus 
expect a modified local law in neighborhoods of this special point. To proceed we need some 
more notation. Recall the domains Si{a,b) in (2.17). For given (small) c ,7 > 0, we set 

S^{a,b) := jz e Si{a,b) : (^jz- 1| > max | a/'1l(ma, Ma), ^JdL{fJ.B, 

and 

5^(a, b) := 5J(a, 6 ) n |z € C : |z - 1| > . (B.4) 

The following proposition presents the local law under the setting (B.l). 

Proposition B.l. Let be as in (B.l), with fixed 4 C o,nd 0. Assume that the 

empirical eigenvalue distributions pA, fJ-B of the sequences of matrices A, B satisfy (2.13). 
Fix any compact nonempty interval X C With the notations and assumptions of 

Theorem 2.5, we have the following conclusions: 

{i) If Pa 7 ^ M/ 3 , then, for any fixed 7 > 0, 

max^|G40-(a7-^B(0)-^h^, 

0a(0 - ^a{z)\ F , |a;^(z) - ujb{z)\ -< , 

hold uniformly for all z G Sx{r]m, 1 )- 
(ii) If Pa = M/ 3 , then, for sufficiently small c > 0 and any fixed 7 > 0, 

max \Gii{z) - {a^ - ujb{z))~^\ F 
l<i<N ' ^ I |z - 1| y/Nfj 

K(0-<^a(0I ^ K(0-wb(0I ^ 

hold uniformly for all z G 1)- 
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Proof. Recall the notation ^ 2 ) in (3.8). In [1] (see the proof of Proposition 7.4 

therein), we proved that under the setting (B.l) and assumption (2.13), one has the fol¬ 
lowing results on the stability of the system (w^, wb, z) = 0: There exists a positive 

constant S such that the following two estimates hold. 

(i): If we have 

{uja,ujb)<S, ze 5 i ( 0 , 1 ), 


and (3.9), (3.14) and (3.15) hold on 5'x(0,1). 
(ii:) If /Tq = we have 

5 


{uja,u!b) < 


2-11 ’ 


2 e 5^(0,1), 


(B.6) 


and (3.14) and (3.15) hold on 5^(0,1), while (3.9) holds on 5^(0,1) with S replaced 



Note that, the proofs of Lemma 5.1, Lemma 6.1 and Lemma 7.1 still work since we have the 
bounds (3.9), (3.14), and (3.15) as well. Although the bound in (3.9) should be replaced by 
in the case /ia = M/Ji it is harmless for our proof. Hence, analogously to the proof of 
Theorem 2.5, one can use Lemma 3.3, Lemma 7.1 and estimates (3.9), (3.14) and (3.15), to 
complete the proof of Proposition B.l. Especially, the proof in the case /ia ^ fip exactly 
agrees with the proof of Theorem 2.5. 

For the case /ia = l^p, we need to replace S by in Lemma 3.3 due to (B.6). In the 

sequel, we simply illustrate the continuity argument in this case. Let z, z' € b), where 

z = E + ir] and z' = E + ir]', with rj' = ij + N~^. In addition, we set zq = z, wi = lua, 
UJ 2 = uiBj ail = and W 2 = in Lemma 3.3. Suppose now that (B.5) holds for z'. Using 
the Lipschitz continuity of the Green function (be., ||G(z) — G(z')|| < A^^|z — z'\) and of 
the subordination functions 01 ^( 2 ) and lob{z) {c.f, (3.9) with S replaced by ■pziy); we can 
choose S in (3.11) to be 


d = 


N'T 1 
|z- 1| 


+ 0(iV-3), 


(B.7) 


In light of the condition k'^ > {c.f., sentence above (3.12), with S replaced by 

p^jj-), one needs to guarantee that SS < \z — l|e, for sufficiently small constant e > 0, 

which is a direct consequence of the assumption that z G S^{a,b) and (B.7). Note that 
||r(z )||2 ^ remains valid since estimate (7.15) does not depend on the stability of the 
system $^^^^e(wA,a;_B,z) = 0, as long as (3.14), (3.15) and (3.9) hold. The remaining parts 
of the proof are analogous to those of Theorem 2.5 and we thus omit the details. □ 
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